Multielectron processes in endohedral atoms at high energies 
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We analyze the role of shake-off and of the final state interactions in inelastic pro- 
cesses in the fullerene shell which follow ionization of the caged atom. We demon- 
strate that in the broad interval of the photon energies the process is dominated by 
the final state interactions. Its contribution is calculated in a model-independent 
way. In the large energy intervals the cross section of the process is close to that of 
photoionization of isolated atom. 

I. INTRODUCTION 

Until now, to the best of our knowledge, the papers on the multiple photoionization 
(mainly the double photoionization) of the endohedral atoms, i.e. of the systems which 
consist of an atom A caged inside the fullerene shell, focused on the case when all photo- 
electrons where ejected from the caged atom . Here we consider another channel in 
which ejection of one electron from the caged atom is accompanied by inelastic multielectron 
processes in the fullerene shell (FS). These may be a single-electron excitation or ionization, 
double ionization, etc. We focus on the case when the energies of the photoelectrons ejected 
from the internal atom E = u — la are large enough 

> 1. (1) 

Here u and are the energies of the photon and of the ionization of the caged atom 
correspondingly (we employ the atomic system of units with e = m = h = 1). We assume 
that the radius of R of the fullerene shell (FS) is much larger than the size of the ionized 
state of the caged atom r^, 

> (2) 
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We consider the spherical fuUerens with the thickness of the shell 

A < i? (3) 

However, the caged atom can be shifted from the center of the sphere. 

This process interferes with the one in which the photon knocks out a FS electron, and 
the latter ionizes the internal atom. This mechanism requires a special direction for the 
momentum of the electron ejected from the FS. The probability is quenched by a small 
factor of the order r\jJ^. Thus we neglect the contribution of this^mechanism. 

As well as in the high energy double photoionization of atoms [3], we can separate three 
mechanisms of the process in endohedrals,i.e. in the systems which consist of atoms A 
surrounded by FS, denoted as A'^Cfq. In the shake-off (SO) mechanism an atomic electron is 
moved to continuum due to its direct interaction with the incoming photon, while the second 
electron is ejected from the FS to continuum by the sudden change of the effective field of 
the atom caged inside FS in the endohedral. In the final state interaction (FSI) mechanism 
the second electron is knocked to continuum due to direct interaction of the ionized atomic 
electron with an electron belonging to FS. In the quasifree mechanism (QFM), which was 
predicted long ago [3] and discovered recently in the experiments absorbtion of the 
photon by the two-electron system takes place almost without participation of the nucleus. 

In the double photoionization of atoms the SO mechanism dominates in the high energy 
nonrelativistic asymptotics. This manifests itself in the behavior of the double-to -single 
cross sections ratio Riyi) = cr^^ {u) / (u) = const. At the energies below the nonrelativistic 
asymptotic both cross section a^^ and its spectrum are the results of interplay between the 
SO and FSI. this is sometimes referred to as the intermediate energy region [5!]. Inclusion of 
the QFM can be viewed as taking into account of the lowest relativistic corrections to the 
photon-electron interaction vertex It leads to increase of the ratio 

We shall analyze the interplay of these mechanisms in the photoionization of the caged 
atom followed by excitations of the FS. Since the QFM requires the coalescence of the two 
bound electrons [6j] , it does not contribute to the considered channel. Thus we consider only 
the interplay between the SO and the FSI. 

Since the FS electrons are separated from the caged atom by the distances of the order 
i? 3> 1, all their interactions and their changes after the photoionization are quenched by a 
factor of the order 1/R, and the probability of the SO is of the order 
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The FSI is determined by the Sommerfeld parameter of interaction between the fast 

n 

electron moving with momentum p and those of the FS [8| 

e = - « 1, (4) 

with p = IpI, c is the speed of hght. Since the binding energies of the valence FS electrons 
Ips < 1, we can write 

e^^« 1. (5) 



We include the FSI terms of the order employing the technique worked out in 10 |. 
Note, however that the probability of FSI contains also the number of the active FS electrons 
N , i.e. the real parameter is .^^A^. In the SO the FS reacts on the change of the field as a 
whole and the probability does not contain the factor A^. 

Because of lack of information about the FS wave function we calculate only the sum 
of cross sections of inelastic processes called also the cross section of absorption, in which 
ionization of the internal atom is followed by ionization of the FS or its transitions to excited 
states, double ionization of FS, etc. 

The cross section of the process with transition of the FS to a particular final state 
n contains the cross section of photoionization of the isolated atom factor. The 

factorization is violated by the terms of the order V/E with V the potential of the FS "felt" 
by the photoelectron. Since \V\ < IRy, such terms can be neglected in our approach. 

o-„ = a^Sn, (6) 

We investigate the behavior of the ratio 

is the cross section of photoionization of the isolated atom in which the photoelectron 
carries the energy E = u — la- We calculate the absorption cross section as the difference 
between the total cross section 

at = (T^St; St = '^Sn (8) 

n 

and the elastic cross section ctq = (J-jSq, i.e. 

r{E) = St{E) - So{E) (9) 
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For the SO mechanism this requires the knowledge of the ground state wave functions 
of the FS. We suggest a simple model which assumes that the FS electrons have a uniform 
distribution inside the shell. The contribution of the FSI can be calculated in a model- 
independent way. It depends only on the number of the active electrons in the FS. We 
calculate the FSI beyond the perturbative approach jo] and trace the energy dependence of 
the function r{E) determined by Eq.Q. 



II. SHAKE OFF 
A. General equation 

We can write for the amplitude of the double photoionization in the considered channel 

F5o = F^($.|^o). (10) 

Here is the amplitude of photoionization of n-th atomic state, \l/o describes ground state 
of the FS electrons moving in the superposition of its self- consistent field and that of the 
internal atom. In the final state $^ one of the FS electrons is moved to the continuum, the 
electrons feel the self-consistent field of the FS and that of the ion with the hole in the n-th 
state of its electronic shell. Here one of the FS electrons is moved to the continuum. The 
matrix element on the RHS of Eq. fllOl) obtains nonzero values only if the angular momenta 
the initial and final states have the same angular momenta. Thus the SO can lead only to 
the monopole transitions. 

The sum of the cross sections of inelastic processes can be written as 

aso = CT,{1 - ($o|^o)'), (11) 

with the cross section of the photoionization, |<l>o) is the ground state of the FS with a 
hole in the electronic shell of the caged atom created by the photon impact. 



B. A model for the ground state 

The wave function \l/o is strongly quenched outside the region 

R<r < R + A, (12) 
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with the radial part depending on 

x = r-R- 0<a;<A. (13) 

The same refers to the function $o- However the values of the parameters R and A in the 
state |\l/o) differ from those in the state 

To estimate the matrix element ($o|^o) in Eg. ( ITT]) we assume that in the ground state 
the FS density does not depend on x. Under this assumptions the wave function of the FS 
electrons with the angular momentum £ is 

^.(r) = {y^'''YU^) (14) 

for r in the interval determined by Eq.( IT3ll . vanishing otherwise {VL is the solid angle). In 
this equation V = AttR'^A is the volume of the FS. 

The FS electrons can be viewed as moving in an effective self consistent field Ui which 
also does not depend on x. The Thomas-Fermi equation 

Po = i^Uof'^„ (15) 
relates the electron density po and the potential 

Uo = li^-'y?^'. (16) 

Here is the total number of active electrons. 

After ejection of the photoelectron a new value of the potential is 

Uf = Uo + Uh, (17) 
where Uh is the potential created by the hole in the state n of the internal atom 

Unir) = [ dh-^p^^ ^ - ^ ^. (18) 
J |r-ra| r R 

where p„ is the electron density in the state n, Ta and r are the coordinates of the atomic 
electron and of that in the FS. 

Thus the final state wave function 

^o{r) = Y,yjyU^) (19) 
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with 

Employing Eqs.( fT4l) and ( IT9l) we obtain 

(^ol^o)-l-^^ (21) 

n 

Assuming, following pjj] the values of the FS parameters for the fuUerene Ceo 

R = 6.02; A = 1.25 (22) 

(in this case = 240 ) we find 

1 - ($o|^o) = 0.046; ($o|^o)^ = 0.91 (23) 
The numerical results do not depend stron gly on the actual values of the parameters. For 



example taking R = 5.75a.u., A = 1.89a.u. 12|] we find 



1 - (<l>o|^o) = 0.060; ($o|^o)^ = 0.88. 

III. FINAL STATE INTERACTIONS IN PERTURBATIVE APPROACH 

A. Lowest order terms 

The amplitude of a process which includes the final state interaction between the fast 
electron and the electronic shell up to the terms of the order is jlO| 

= F^'^ + Fi'^ + F^^\ (24) 

where the upper index denotes the number of interactions between the fast electron and the 
FS, index x labels the final state of the FS. The amplitudes fJ*^ contain the amplitude of 
ionization of the isolated atom factor Q 

FP = F^T^^; 2 = 0,1,2. (25) 

Here T'^°) = is the SO matrix element-see Eq. ffTOj) . The accuracy of this equation is 

Va/E with Va the potential energy of the photoelectron in the field of the caged atom. 
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One can write for the cross section of the process with transition of the FS to a particular 
final state x [lO|(see Eq.([6])) 

= iTjof + 2rio)i?eTii) + l/mTi^f + 2T^°^ReTi^\ (26) 

Now we calculate the FSI amplitudes. One can write 

TW = (<|..|f/i|vI/o), (27) 

where Ui = ?7i(r*^'^^), with k labeling the FS electron, t/i(r*^^)) is its interaction with the 
photoelectron in the lowest order of the FSI. One can write 

t/,(r(^)) = 1 1 ^G(f),(/)e^(^r<^'), (28) 

where G{f) = 2(/i^ — (p + f)^ + it^) with fi"^ = p"^ + 2efi is the free electron propagator, 
energy Efi is transferred by the FSI, g{f) = Ati/{P + A^), A — )■ 0. Keeping only the term 
proportional to the large momentum p in denominator of the electron propagator we can 
put 

G(f) = ~^ . , (29) 
2(p ■ i) — IV 

Using the well known formula 

11/"^ dx 



a c 



(^ax + c(l — x) j 



(30) 



we obtain 

dx 



2(p-f)-iz/ /2 + A2 Jo (2{p ■ - x) + Px + X^x - iu)^' 
Introducing y = [1 — x)/x, f = i + py and integrating over f by using the relation 



(31) 



d^f 47re*(f-'') 1 d e 



ibr ^ibr 



(27r)3 (/2_52_^^)2 2hdh r 2b 



b'=pV-X' (32) 



we find 
this leads to 

Ui=i^J2^n{r^^^ ~rf^)\ (34) 
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Thus the amphtude T*^^^ is mostly imaginary 

r« = t^{^.\J2Hr^'^ - )A|^o). (35) 

k 

The imaginary part dominates since the pole of the electron propagator provides the leading 
contribution. This means that the photoelectron passes the distances of the order of the 
FS radius R ^ 1, and interacts with the FS electrons at the region of their location. The 
divergence at A = is just the Coulomb phase of the interaction between the photoelectron 
and the electronic shell jlO|. The divergent contributions will cancel after the second order 
terms will be taken into account. 

Since the leading contribution to Tx^^ is imaginary, while Tx^^ is real, the leading non- 
vanishing contribution of the FSI is of the order In order to find it, one has to include 
the second order amplitude tJ^"* and the terms of the order //p in the first order amplitude 

7^(1) 

-t X ■ 

The second order amplitude can be written as 

Ti^) = ($.|?72|vl>o), (36) 



with 



^2 = ^ E / + f2)^(/2)e^(^^-'^^')e^(^^-'^^^ (37) 

Using Eq. (j29ll for the Green function G and putting in the integrand 

1 1 1/1 1 1 1 \ , , 



p • fi) (p • (fi + f2)) 2 V(p . fi) (p . (fi + f,)) (p . f,) (p . (fi + f^)) 



2(p.fi) (p. fa)- 



we find that 



U2 = Ul/2 (39) 



Here again the FSI takes place in the region of the FS. As to the terms of the order f /p in 
the amplitude they determine its the real part ReT^\ they describe the interactions 
between the photoelectron while it is close to the caged atom and the FS. Thus they have 
additional factor 1/R <^ 1. Thus the second term on the RHS of Eq.f l26|) is 

Tfi?eT«-^|($.|vI^o)r-^ (40) 
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Hence,neglecting the terms of the order in the FSI terms, we find St = 1 while 

So = |Tf p + |(vI/o|f/i|^o)P - (^olf/i'l^o). (41) 

We can put r^'^^ = R in Eq. (lM|) . As expected, the terms containing A cancel on the RHS of 
Eq. (l4Tll . Hence we can employ 



f/i = i^A; A = ^ ln(l - t^'^)) (42) 

k 

with t'^^^ = rf^ jr^^). Thus Eq. (l4Tl) takes the form 

^0 = |ro(°^r + e'l(((^o|A|vl>o)P - (v&olA^lv&o)). (43) 

Here the three terms on the RHS describe the SO, FSI and their interference correspondingly. 
Direct calculation provides 

r{E)=C + eN; C = l-|To^°^r (44) 

Employing of closure requires that the energy E is large enough to include all important 
excited states 

^ > £. (45) 

At large energies e ^ Ips the energy distributions drop as and thus the energy losses 
e are determined by Ips ^ e ^ E. They are [13 1 

(^N ^ E 

Some additional data can be obtained by studying the distributions of the electrons, 
ejected from the FS. 

B. Partial wave analysis 

The ratio r{E) given by Eq. (l44l) can be written as the sum of the contributions of the 
partial waves 

nt n£ n£ 

with the sum carried out over the excited states n 7^ 0. 
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The first and the third terms term on the RHS obtain nonzero values only for i = 0. 
Expanding in terms of the Legandre polynomials 

ln(l-t) = ^a,P,(t); a, = / dtPi{t)ln{l - 1), (48) 

we see that the contributions of the terms with £ = to the third and second terms cancel. 
Thus 

r{E) = 5^(5,0 + eN biAni)', (49) 



i=i 



with 



while 



Or, employing closure 



Presenting 



iE) = l-Al, + eNY,h. (52) 



2i + l 



(53) 



£2(£+l)2 £2 (£+1)2 

we see that indeed ^ 6^ = 1. 

The dipole term which dominates in the FSI provides 3/4 of its contribution. The am- 
plitude of the dipole transition to a particular excited state n is proportional to the overlap 
of the radial wave functions ($]!^]^|\1/q). Thus investigation of the spectrum of the electrons, 
ejected from the FS would provide the data which is cornpliment ar y to that obtained from 
the studies of direct photoionization of the fuUerene Cqq |14|.|15|. 

Note, however that the perturbative results, presented in the last two Subsections are 
true if ^^A^ ^1. In the case of the fullerene Cqq this means that E ^ 5 keV. At such 
energies all = 360 electrons of the FS can participate on the process. In order to analyze 
the lower energies we must go beyond the perturbative approach. 



IV. FINAL STATE INTERACTIONS BEYOND THE PERTURBATIVE 

APPROACH 
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One can see that Eg. (1381) can be generalized for the case of arbitrary number n of inter- 
actions between the photoelectron and the FS. Introducing a„ = (p • (fi + f2 + ..f„)) we can 
write 

1 . 1... . i = li. (54) 

Oi a2 an n\ a" 

This equation which can be proved by the induction method, was used for calculation of 



161]. Thus the amplitude of 



the radiative corrections in the electromagnetic interactions 
transition to the state {(^x\ is 

= ($,|e*«^|^o), (55) 

with A defined by Eq.g2]). Thus 

t{E) = 1 - |($o|e^«^|^o)|' ^ 1 - |(^o|e^^^|vl/o)r (56) 

The last equality is due to the fact that the change of the field of the caged atom provides 
the perturbation of the order 1/R. Hence 

r(E) = l-|(<l>o|n,(l-t('=))^«|v^o)r- (57) 

This provides 

r{E) = 1 - ^^:^l($o|^o)r = 1 - e-'^'-^'^^^^l^o)]'. (58) 
Thus we can write 

r{E) = C + rpiE); r^(E) = 1 - e-^i-(i+€'). (59) 

Here C defined by Eq. fH^ is the SO contribution which does not depend on E, rpi^E) is the 
contribution of FSI. 

If the photon energy is so large that A^^"^ ^ 1 ( i.e. E ^ 300 eV) we find 

rpiE) = 1 - e-^^\ (60) 

At these energies A^ is the total number of the FS electrons, i.e. A^ = 360. 

V. TOTAL CROSS SECTION 

Now we can trace the energy dependence of the cross section of absorption by the FS in 
this process. Employing Eqs. (H5|l and (H6|) we find that closure can be used at > 50 eV. If 
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the energy E is smaller than the ionization potential of the core Is electrons Ic ~ 315 eV., 
the ratio rp{E) is determined by Eq. fl59|) with N = = 240 the number of the valence 
electrons. At E close to Ic we find 1 — rpiE) ^ 2 ■ 10~^. 

At E > Ic the core electrons are involved into the process as well. While E is of the 
order Ig their contribution can not be calculated by employing closure since some of the 
exited states can not be reached due to the energy conservation low. However, Eq. (159!) with 
N = = 240 provides the lower limit for the value of rp{E) at these energies. 

At larger energies E ^ Ic the ratio rp{E) is determined by Eq. fl60l) with the core electrons 
included, i.e. = 360. ki E = 2 keV we find 1 - rF{E) ^ 0.09, i.e. rF{E) is still very 
close to unity. At = 5 keV rp^E) ^ 0.62, dropping as 1/E at larger energies, following 
Eq.( l44l) . The FSI and the SO contributions to the ratio r{E) Eq.( l59i) become of the same 
order at ii^ > 50 keV. 

If the binding energy of the ionized state of the internal atom and the photon energy 
are small enough, both SO and FSI contributions to the cross section are enhanced by the 
same factor. This happens because external photon is strongly infiuenced by the FS due to 
polarization of the latter. This effect manifests itself in a factor 

D{u) = 1 - a{uj)/R^ (61) 



in the amplitude of ionization of internal atom [17|. Here a{uj) is the dynamic dipole 
polarizability of the FS. Since a{u) < the factor |D(a;)p increases the cross section of the 



single photoionization 17| 



<oK^) = ^H^)D\u). (62) 

The polarizability a{uj) has a strong maximum due to FS giant resonance at w ~ la.u. 
The role of polarization diminishes with increasing of the photon energy. It becomes 
negligible at Wmax ~ 2.5a.u, thus here D'^iu) ~ 1. However, the characteristic binding energy 
of a FS electron is / = 7eV. If the binding energy of the ionized state of the caged atom 
is of the same order or smaller, the factor -D^(w) increases both SO and FSI contributions 
to the double photoionization cross section at the lower limit. Note that the factor -D^(w) 
enters the cross sections of the double and single photoionization in the same way. Thus it 
cancels in their ratio r(i?)-Eq.(I71). 
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VI. SUMMARY 

We investigated the high energy photoionization of the endohedral atom A@Cqq followed 
by inelastic processes in the fuUerene shell. We traced the energy dependence of the ratio 
r{E) of the cross section to that of photoionization of isolated atom. In a broad interval of the 
values of the photoelectron energies E the ratio is dominated by the final state interactions 
and is calculated in a model independent way beyond the perturbative approach. The ratio 
is shown to be close to unity until we reach the region of E of several keV. At the energies 
higher than E = 5keV where r{E) = 0.62 it drops as 1/E. The contributions of the FSI 
and of the shake-off become of the same order at very high energies E > 50 keV. 

The work was supported by the MSTI-RFBR grant 11-02-92484. One of us (EGD) thanks 
for hospitality during the visit to the Hebrew University of Jerusalem. 
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